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A mean field calculation of the T — > limit of the Hall conductance of electron-doped cuprates 
such as Pr2-xCexCu04,+s is presented. The data are found to be qualitatively consistent with 
the reconstruction of the Fermi surface expected upon density wave ordering. The magnitude of 
the density wave gap is found to be large. The Hall resistance exhibits a nonanalyticity at the 
quantum critical point for density wave ordering, but the amplitude of the anomaly is found to 
be unobservably small. The quantum critical contribution to Rh{B) is determined. Quantitative 
discrepancies between calculation and data remain, suggesting that the experimental doping is not 
identical to the Ce concentration x. 
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I. INTRODUCTION 

The unusual properties of high temperature cuprate 
superconductors continue to challenge our understanding 
of the physics of electrons in metals. The related ideas 
of "quantum criticality" ^ and density wave (or "stripe" ) 
order and fluctuationsSii have been much debated, but 
in the hole-doped cuprates the discussion has been in- 
conclusive. Further, while it has been proposedii^ that 
many of the anomalous properties of the materials could 
be explained by a quantum critical point at or near "op- 
timal doping" , there is no unambiguous evidence of long 
ranged order setting in near optimal doping in the hole- 
doped cuprates. 

The electron- doped cuprates offer a new perspective 
on these issues. In these materials commensurate (tt, tt) 
spin density wave order has been detected by muon spin 
rotation^ and neutron scattering measurements The 
order exists over a wide range of dopings and vanishes at 
a critical doping Xc near the "optimal" doping x = 0.16. 
Thus, in this material an unambiguous quantum critical 
point exists, separating an apparently disordered phase 
from a phase with a well-defined ((tTjTt) commensurate 
magnetic) order. It is therefore of interest to examine the 
effect of the order and criticality on material properties. 
In this paper we consider the low temperature Hall resis- 
tance, which is sensitive to the rearrangment of the Fermi 
surface caused by the onset of the spin density wave or- 
der. We apply to the electron-doped cuprate case a mean- 
field analysis presented previously in the Iiteraturej2i2ii£ 
and compare the results to experiments. ^^'^^ 

The rest of this paper is organized as follows. Sec. UTI 
defines the model to be studied. The formal solution to 
Boltzmann equation with relaxation time approximation 
for 2-dimensional systems is derived in Sec. IIIII Sec. IIVI 
gives the numerical results for Rh in weak magnetic field 
and compares them to experimental data. In Sec. Ivl we 
analyze our formal solution in the critical region, x < 
Xc- The sensitivity of our numerical results to variation 
of parameters is studied in Sec. IVII The possibility of 
explaining the x > Xc experimental results by calculating 
the real part of the self energy due to spin fluctuations 



is investigated in Sec. IVIII We conclude our calculations 
and discuss our results in Sec. IVIIII 

II. MODEL 

We consider electrons moving in a square lattice of unit 
lattice constant with dispersion 

Sp — —2ti (cOSPx -\- COSPy) -t- 4^2 cos Pa; COSpy 

— 2^3 (cos 2px -f- COS 2py) (1) 

The "canonical" values of the band parameters are ti = 
0.38eV, t2 = 0.32ti and = 0.5*2- As will be shown 
below, our main results are not sensitive to the precise 
parameter choices. We are interested in electron dopings, 
corresponding to two dimensional carrier densities per 
unit cell n = 1 + X > 1. 
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FIG. 1: Fermi surfaces calculated from Eq \^ for x — 0.15, 
band parameters ti = 0.38eV, t2 — 0.32ti, = 0.5*2 and 
backscattering ("gap") values A = (dashed line), A = 0.2eV 
(light solid line) and A = OAeV (heavy solid line). The dotted 
line is given by Sp+Q — /i(A = 0). 

We assume that for x less than a critical value 
commensurate spin density wave (SDW) order occurs, so 
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the electrons are subject to a coherent backscattering of 
amplitude A and wavevector Q = (tt, tt) which doubles 
the unit cell, changes the dispersion to 
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sp + ep+Q ± V (£p - £p+Q) + 4A2 (2) 



and reconstructs the Fermi surface. If A = then (at 
the relevant dopings x) the Fermi surface consists of one 
large hole surface centered at the X point (tt.tt) (dashed 
line in Fig ^) . If A 7^ then the Fermi surface is re- 
constructed. The details depend on the band filling and 
the magnitude of A. If the chemical potential is such 
that A < A*(/i) = £(7r/2,7r/2) — /X = 4^3 — then both 
bands cross the Fermi level and the Fermi surface in- 
volves two symmetry-inequivalent hole pockets centered 
at (7r/2, ±7r/2) and one electron pocket centered at (0, tt). 
One of the hole pockets and two portions of the electron 
pocket are shown as the light solid lines in Fig^ How- 
ever, if A > A*(/i), the lower band is completely filled 
and only the electron pocket remains (heavy solid line in 
Fig^. For the parameters used in FigQ] A* « 0.26eV. 

As X is decreased below Xc we assume that the gap 
magnitude increases. However, the band filling and hence 
the chemical potential decreases. Therefore, depending 
on the relative growth (with decreasing x) of ^ and A 
either the hole pockets vanish at a critical a; = xi, or it 
persists to a; = 0. The comparison to data given below 
suggest that the former circumstance obtains in the ac- 
tual materials, and that indeed the gap opens so rapidly 
that xi is of order 0.12, not much less than Xc ~ 0.16. 



III. CALCULATION OF HALL RESISTANCE: 
FORMALISM 

We require the T transport coefficients. In the 
T ^ limit the only scattering is due to impurities, 
and we assume a Boltzmann equation description is ad- 
equate, so that the physics is described by a momentum 
space distribution function g(p) = f{p) + hp with f{p) the 
Fermi-Dirac distribution. We also adopt the relaxation 
time approximation, linearize in the applied electric field 
E and take the T ^ limit so that for each band b, the 
Boltzmann equation becomes (e = |e|): 



eE ■ v'SiE; - fi) 



d_ 

dp' 



(3) 



We are interested in two dimensional electrons. We take 
B \\ z and note that B-{'if^x -^) corresponds to a deriva- 
tive along the Fermi contour. Defining s to be the coor- 
dinate along the Fermi arc, noting that the structure of 
the equation implies 



(4) 



and temporarily dropping, for ease of writing, the band 
superscript b allow us to rewrite Eq [3] as 



with (I = VT is the mean free path, and $o = hc/2e is 
the superconducting flux quantum): 



a{s) 
Eis) 



-{s)B\v{s)\ _ Trl{s)B 



he $0 

eT{s)E ■ V 



(6) 
(7) 



Eq [S] may be formally solvedi^ in terms of the Green 
function K{s, s') satisfying 

[l/a{s) + d,]K{s,s') = S{s-s') (8) 

with boundary condition 

K{s^ s'+,s')-K{s^ s'-,s') = l (9) 



h{s) = (b ds'K{s,s')E{s')/a{s') 



(10) 



The currents longitudinal and transverse to E may then 
be constructed by averaging over the velocity as usual. 
Assuming a square lattice, restoring the factors of veloc- 
ity and the band label (if necessary), and symmetrizing 
or antisymmetrizing for axx,xy, we find 



dsds' 

^ Pb 

dsds' 



K{s,s')cos[(f,{s) - (l){s')] (11) 



CTxy = ctq * — —K{s, s') sin[0(s) - 0(s')] (12) 
J Pb 

where ctq = e'^/h is conductance quantum, pB = ^tt/Ib 
is proportional to the inverse of the magnetic length I b = 
(<i>o/7ri3)^/^, and (j) is the polar angle of Fermi velocity 
on the corresponding Fermi contours. The formulae are 
valid as long as magnetic breakdown can be neglected. 
A criterion for this is given below. 

The problem has several momentum scales. Obtain- 
ing a non-vanishing Hall coefficient requires breaking of 
particle- hole symmetry, either by a variation of scattering 
rate, a Fermi surface curvature, or a variation of velocity. 
All of these phenomena are characterized by momentum 
scales of the order of the size of the Brillouin zone; we 
refer to these scales generically as pq. We take these to be 
large, so we may neglect Landau quantization. A second 
important scale pA = A/vp is set by the SDW backscat- 
tering amplitude. A, and vanishes as A — > 0. A third 
scale is the range, a, of the kernel, K (Eq [Hj). We dis- 
tinguish weak field (a < pa) and strong field (a > pa) 
regimes. A fourth significant scale is the inverse of the 
mean free path, l^^. The results we present are valid 
for Pq ^ PA,a,l~^. In addition, we expect disorder to 
modify the SDW behavior significantly for paI ^ 1; so 
our results do not apply for A < vp/l. In this case we 
may neglect the periodicity and write 



[l + a{s)ds]h{s) ^E{s) 



(5) 



K{s, s') — exp 



dx 
a{x) 



e(s-s') 



(13) 
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FIG. 2: Local view of the Fermi surface in the vicinity of one 
of crossing points, Sp = Ep+Q = /i, shown in Fig p, sp 
are respectively the Fermi Imes for the normal state and the 
normal state Fermi line shifted by Q. The curvature has been 
neglected. + denotes the electron pocket of SDW state and — 
denotes the hole pocket. Arrows along Fermi lines show the 
circulation direction determined by B x v with B |[ — z. They 
also indicate the direction of increase in arclength coordinate, 
s. V , V'' , and V are Fermi velocities along corre- 
sponding Fermi lines. The dashed line is the diagonal of the 
first quadrant of Brillouin zone connecting (0, tt) and (7r,0). 
The dot-dashed line is perpendicular to the dashed line: if 
the normal state energy dispersion is linearized, and the mo- 
mentum dependence of velocity is neglected, the problem ts 
symmetric under refiectton through the dot-dashed line. 

Finally, we note that all of the formulae are derived on 
the assumption that a state evolves along a given Fermi 
contour and does not jump to another one, i.e., that mag- 
netic breakdown may be neglected. By applying the con- 
siderations of Ref. [Tj to the model studied here, we find 
the condition for absence of magnetic breakdown is, in or- 
der of magnitude, that the inverse of the magnetic length 
Pb — (7ri3/$o)^^/^27r be less than the gap momentum 
scale PA- Our results are only valid if a,pA > Pb- The 
condition that breakdown occurs at a field higher than 
the weak-strong field crossover is p^l > 1, which as noted 
above is a necessary condition for the validity of the ap- 
proach. 

IV. WEAK FIELD LIMIT: EVALUATION AND 
COMPARISON TO DATA 

In this section we analyse in detail the case in which 
the field is weak enough that a is small relative to all 
scales in the problem. In this case the s' integrals in Eqs 
|llll2j are dominated by s' ~ s and may be performed, 
leaving 

<Jxx = '^Q^^ j' dsl{s) (14) 
^.y = OQ^jds{l{s)/lBf^ (15) 

for the longitudinal and Hall conductance per plane. Eq 
[Tlj is exactly the Drude result (average of mean free 



path over the Fermi arc), while Eq ^B] reproduces Ong's 
resuliii^ 

B \ f -> -> 

CTxy = CTq — — (b dl X I - z/2 (16) 

Further analysis requires information about the scat- 
tering rates. The two commonly made assumptions are: 

• mean free path is constant over the Fermi surface 

• scattering rate is constant over Fermi surface 

The constant scattering rate assumption may be derived 
theoretically by assuming weak (Born approximation) 
point-like scatterers. Constant mean free path is more 
appropriate for very strong (unitary) scattering. 

Ongi^ shows that the constant mean free path approx- 
imation implies that the Hall conductance for each Fermi 
surface pocket is a constant, depending only on the value 
of the mean free path and the sign (electron- or hole-like) 
of the Fermi pocket. Therefore, in the constant mean 
free path approximation, axy does not change as x is de- 
creased through Xc- The reconstruction of the Fermi sur- 
face means that one goes from one hole pocket {cTxy oc P) 
to two hole pockets {a^y oc 2P) and one electron pocket 
{axy oc —P) so that the total cjxy (x P{2 — 1) = P. 
However, at x = xi the hole pockets vanish, and in the 
constant mean free path approximation the Hall coeffi- 
cient would jump discontinuously from a value close to 
the Rh ^ 1/(1 ~ 2;) implied by the initial, large, hole- 
like surface, to the value Rh ^ —1/x implied by a single 
electron-like pocket. This behavior is in evident disagree- 
ment with data,^^'^^ which indicates instead a smooth 
change beginning at x ~ 0.16, which is approximately 
the doping at which antiferromagnetism sets ini 
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FIG. 3: Solid line: Hall coefficient calculated from Eqs 
for parameters ti = 0.38eV, t2 = 0.32ti, t2 = 0.5t2 and 
l^{x)[eV] — 0.6y'l — x/0.16. Long dashed line: Hall coef- 
ficient calculated for parameters ti = 0.38eV, t^ = 0.32ti, 
t2 = 0.075eF and A{x)[eV] = 0.7^1 - 2:/0. 16. They are 
calculated with constant 1/r approximation. Dotted line: 
curve Rh = —1/x for x < 0.16 and Rh = 1/(1 — x) for 
X > 0.16. Stars: experimental datam^ Solid squares: exper- 
imental datar— The open box shows the approximate size of 
the region for which the weak field-strong field crossover may 
be observable with present accessible fields. 
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We therefore turn to the constant 1 /t approximation. 
In this case, a^y for each pocket is weighted by the square 
of the Fermi velocity for that pocket and more varia- 
tion occurs. The solid line and long dashed line in Fig 
Oshow typical results obtained by evaluating Eqs |14ll6j 
in the constant scattering rate approximation and assum- 
ing a backscattering amplitude growing proportional to 
y^Xc — X. One sees that for x > Xc the Hall resistance is 
hole-like and depends only very weakly on x. At x = Xc 
a nonanalyticity occurs, and Rh begins to decrease. For 
the parameters chosen, xi ~ 0.12. At this doping the 
hole pockets vanish, a further change in curvature occurs, 
and Rh approximates very closely the 1/x behavior ex- 
pected from the single remaining electron pocket (shown 
as the dotted line). 

FigElalso shows, as stars and heavy squares, the exper- 
imental resultsiii*i^ These data are in qualitative agree- 
ment with calculation, but important quantitative differ- 
ences occur. At the lowest doping (nominal carrier con- 
centration X = 0.11 for starsii and x = 0.05 for heavy 
squareSiiS.) the measured Rh is substantially below both 
the calculated value and the 1/x curve. It seems to us es- 
sentially impossible to obtain a Hall resistance satisfying 
\Rh\ > 1 1/x I; we therefore suggest that the nominal dop- 
ing level, X, differs from the actual doping. The presently 
available data do not allow us to locate the point xi with 
precision, but strongly suggest that xi > 0.1. 



V. CRITICAL BEHAVIOR 

In this section, we study in more detail the behavior 
of Hall conductance axy in the vicinity of a; = Xc- When 
the system is tuned through the critical point (which we 
believe to be Xc ~ 0.16 in the electron-doped cuprates), 
the coherent backscattering introduces reconstruction of 
Fermi surface, changes its local curvature, and thus af- 
fects the Hall conductance. The weak field limit was pre- 
viously studied by Bazaliy et in the context of the 
3d SDW compound Cr. They obtained a result equiva- 
lent to our Eq |3H and a qualitative result equivalent to 
the weak field limit of our Eq |22] . And the strong field 
limit was studied by Fenton and Schofield^ who found 
a result equivalent to our Eqs |32I33| . The two dimen- 
sional nature of the present problem allows us to carry 
the analysis further, obtaining prefactors for axy in the 
weak field limit and a treatment of the weak to strong 
field crossover. 

To study the vicinity of the critical point, we use Eq 
[T^ to write an equation for the change in the Hall con- 
ductance as the gap is opened, Sa^y = <7xyi^) — o'xyi^ = 
0). For small A the important changes in the Fermi sur- 
face occur near the points where the A = Fermi surface 
crosses its translation by Q, i.e. where Sp — £p+Q — /i. 
One of these crossing points is shown in Fig [21 The 
change of chemical potential is second order effect of A, 
and thus neglected in our analysis. To simplify notation, 
in this section, we measure /i from its value at critical- 



ity. Thus the SDW Fermi surface with backscattering 
potential A, = 0, can be written as SpSp+Q = A^. 
Any change in the scattering rate is similarly of order 
A^ and is neglected. Near such a crossing point we may 
write Saxy in terms of the contributions from the two 
SDW bands and the original Fermi surface, p, and 
its image under translation by Q, sp, as 

Sdxy ^ / dsds' i ^ Kb{s, s') sm[(j>i,{s) - (j>b{s')] 

a—p.sp ^ 

We shall see that this integral is dominated by s,s' close 
to the crossing point, with "close" determined by the 
larger of a (evaluated at the crossing point) and pA, so 
that we have neglected the fact that the Fermi surface 
closes on itself and have extended the integration to in- 
finity. 

We analyse Eq jl7| by expanding the energy disper- 
sions about the crossing point. In the first approxima- 
tion, one may consider taking r to be constant and lin- 
earizing the dispersions about the crossing point, writing 
Ep — ■ Sp and Sp+q = ■ 6p, with v^'"'''' constant. 
In this approximation, the A = Fermi surface and its 
image are straight lines, so are independent of s, 

so the A = contribution vanishes, and the two SDW 
bands are mirror images of each other with respect to 
the dot-dashed line in Fig |21 so K+{s) = K-{s) and 
4>+{s) — 0+(s') = — (0_(s) — 0_(s')), so that Saxy van- 
ishes. Obtaining a non-zero Saxy requires consideration 
of effects which break this mirror symmetry. These in- 
clude variation of t or with momentum as well as 
curvature of the A = Fermi surface. Formally, all of 
these effects are of order p^^ , and we work to leading 
order in this quantity. Then we may write 

a^{.s) = a{s)[l±Sai.s)] (18) 

where a(s) is the corresponding function in presence of 
mirror symmetry, and Sa(.s) represents the broken mirror 
symmetry, and thus is of order O{l/po). In terms of these 
functions: 

K^{s,s')^K{s,s')(l± rdx^-^) 

V Js' a(x) J (19) 

^ K{s,s'){l±SK) 

We also have: 

0+(s) - cj,+ {s') + 0_(s) - 0_(s') = f{s,s') (20) 

where / is also of order l/po. Thus, Eq ^21 can be ex- 
panded as: 

ScTxy ^ J dsds'^2K{s,s')f{s,s')cos[(t)+{s) - (?!)+(s')] 

4- 2K{.s, s')SK{s, s') sin[(/)+(s) - 0+(s')] 

- Kp{.s, s')ff(s, s') ~ Ksp{s, s')g{s, s')| 

(21) 
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with g{s,s') of order 1/po- In obtaining this equation, 
we have noticed that /(s, s'), SK and g{s, s') are aheady 
of order 1/po: so that aU the remaining quantities can be 
taken from their unperturbed counterparts. We also no- 
tice that a^(s) — a'^P(s), and in this approximation they 
are constant, denoted by in the following discussion. 

The detailed analysis of this equation is somewhat in- 
vovled, except for the strong magnetic field limit to be 
discussed below, partly because it is not easy to obtain an 
analytic relationship between the arclength coordinates 
s and momentum coordinate-p for SDW Fermi line. We 
shall argue that this Sa^y has a scaling form: 



(22) 



such that -F'(O) = 1 and F{x oo) cx x, with the prefac- 
tor C of order 1/poj which will be computed in the special 
case of constant relaxation time in the following, and the 
exact form of F depending on the manner in which the 
symmetry between + and — branches is broken. 

It has been noted by Bazaliy et a.'^^ that in the vicinity 
of the crossing points, p can be parametrized by £p and 
£)3_i_q: first expanding the dispersions 



Ep = vP* ■ Sp + -mij6p,Spj + 0{5p^) 



v'P* • (5p + -n,j6p,Spj + 0{Sp^) 



(23) 



where, dp ^ p ~ p\ wP^'p* - vP'^P(p*) = (VpE^.p+g)*, 
niij = {^'^ep/^p^^pj)*, = {d'^Ep+q /dpidpj)* , and the 
superscript * stands for the corresponding value evalu- 
ated at the crossing point, and then inverting them 



6p = UiEp + U2£p 



A2 

(ui + -^U2)ep (24) 



with Ui 
From Eq 



v°P*x[vP*xv^P* 



p^yi and U2 



yy. X [V°'" XV^" 

(vP* xv=P*)^ 



dp = uidsp + \i2dep+Q = (ui -U2)dep (25) 



where the second equality of Eqs |24I25| only holds on 
the SDW Fermi surface. Thus, 



ds 



dsr, 



I [(v'Pf + (vPf ( —V 

+ (26) 



a; 

£2 
- 

And the SDW Fermi velocities are: 



v±(p) 



A2 .2 



rv^(p) 



P .r^P 



£2 + A^ 



v^P(p) (27) 



Now return to Eq [22 . a(s) = TrSrof (s)/$o, where tq 
is the relaxation time at crossing point, v{s) is the SDW 



Fermi velocity when mirror symmetry is present, in which 
case vP and in Eqs |2()l27j are constant vectors. Sub- 
stituting this and Eq |2S1 into Eq |23 , scaling energy vari- 
ables in unit of A, and defining B — B{vpTTTo sin A(/)/$o), 
where vp = vP'^p and A0 is the angle made by the nor- 
mal state Fermi velocities, and v''^', shown in Fig[21 we 
can obtain the above scaling form Eq |221 • In the follow- 
ing, we shall discuss Eq |22j in two limits: weak magnetic 
field, aP <C PA, and strong magnetic field, aP ^ p^. 

In weak magnetic field, we can start from Eq ^^l- Both 
vP and r can vary in momentum space. To lowest order, 
these effects are additive. If we keep the p-dependence 
of relaxation time r, assuming its dependence has form: 



Tp ^tq[1 + A{ep + Ep+q)] 



(28) 



and keep velocity of normal state, and v^p, to be 
constant, Eq can be evaluated for each Fermi pocket, 



givmg 



da^y = aq4:TQA——z 
5>o 



(29) 



In the case in which the normal state velocities and 
v^P are p-dependent and relaxation time t is constant, 



differentiating Eq 
defining vectors: 



substituting Eq [21] into it, and 



f]P = (miiwij; -I- mi2Uiy)± + {m2iuix + m22Uiy)y 

J% = {miiU2x + mi2U2y)x + {m2lU2x + m22U2y)y 

= {niiuix + ni2Uiy)x + {n2iUix + n22Uiy)y 

rff' = {niiU2x + ni2U2y)x + {n2lU2x + "22W2y)y 

we then have v'' = v''* + ^v'' = v*** + {rfl + rj^{A/sp)^)ep 
and dv'' = (yyf - '4{A/sp)^)dep, {b = p, sp). Eq ^ then 
gives: 

R A 

Saxy = CTqT'^z-M + 7tP + 3rg + 3fJ^P) x (v^^*- v^^)] 

(30) 

confirm 



These weak magnetic field results, Eqs 
our scaling form for Sffxy, Eq |22j . in the limit aP /p^ <ti 1. 
Also the prefactor C can be calculated from these two 
equations. It would be of interest to study Eq |3Uj in more 
detail, since all the parameters in Saxy/^xyix — Xc^-) 
are set by normal state energy function Eq Taking 
Cxyix = Xc+) from numerical calculation, daxy/cTxy{x = 
Xc+) = — 6.1A, shown as solid line in the upper left panel 
of Fig ^ which is compared with the results of numeri- 
cal evaluation of Eq ^B] (solid circles in this panel) . In 
weak magnetic field, the longitudinal conductivity axx 
has been treated in great detail for the case of the 3d 
SDW transition in Cri^ Following essentially the same 
procedure, we obtain a similar expression for 6axx- 



T (vP* 
ocTxx — —crq — 



A 



(31) 



Taking axx{x — Xc+) from numerical calculation, 
S'^xx/o'xxix = Xc+) = — 2.1A, shown as solid line in up- 
per right panel of Fig^ compared with the results taken 
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FIG. 4: Upper left: numerical results for (Jxy{^) /<Jxy[^ = 0) 
as a function of A near critical doping (solid circle) and func- 
tion 1.0 — 6.1A (solid line). Upper right: numerical results 
for (Ta;2,(A)/(Ta;2,(A = 0) OS function of A (solid square) and 
function 1.0 — 2.1A (solid line). Lower left: numerical re- 
sults for Rh as a function of A (solid diamond) and func- 
tion — 2.4A (solid line). Lower right: numerical results for 
Rh as a function of x near critical doping. The param- 
eters are: ti = 0.38eV, t2 = 0.32ti, t-i = O.OTSeV, and 
A{x)[eV] = 0.6^1 - x/0.16. 

from numerical evaluation of Eq |14) (solid squares). The 
numerical results of Rh = cTxy/iBa"^^) are shown as 
solid diamonds in the lower panels of Fig 01 The re- 
connection of Fermi surface suppresses both longitudi- 
nal and transverse conductivity. In the vicinity of crit- 
icality, Rh/Rh{x = Xc+) « 1 -I- Saxy/uxyix = Xc+) - 
'^So'xx/o'xxix = Xc+) ~ 1.0 — 1.9A. Since A ~ y/xc — x, 
we expect a square root behavior in R^ (x) . There actu- 
ally is, as shown in the lower right plot of Fig^J however, 
this regime is so small that it may be difficult to observe 
in practice. 

In strong magnetic field regime a^/pA S> 1, the inte- 
gration Eq |2J is dominated by (s > 0, s' < 0), with 
f{s,s'), g{s,s') {s - s')/pq, Sa{s) sgns • s/po, and 
(j)+{s) - (j>+{s') w Acj), then 

Saxy = aq ^^^^ [cos A0 - 1 + a sin A0] (32) 
tt^IbPo 

with a ~ ^ 1 and / the mean free path in normal 

state. This result, Saxy ^ and independent of A, con- 
firms our scaling function for aP jpt^ » 1. To justify this 
result, let's consider a simple case in which the normal 
state velocities v^'''*p are constant, while relaxation time 
is assumed to be of form: Eq |2H1- Under this condition, 
the only contribution to 8oxy in Eq comes from the 
second line containing 8K. Changing to variable, it 
gives 

Scyxy=cjQ'-^BAf{BI^) (33) 

]{x) is a function of only one variable, numerical eval- 
uation of which gives the solid line in Fig [Sj along 
which shown as dashed line is the linear function j(x) = 
2a;, representing the calculated asymptotic behavior as 
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FIG. 5: The scaling function /(B/A) for Saxy under the lin- 
earized dispersion, constant velocity approximation. The unit 
of B/A IS B/A. 



X ^ oo. The prefactor 2 is calculated by relating 
a in Eq 1221 and ^ in Eq |2H1 in the hmit A ^ 0: 
a/po — A{vf sin Ac/)) /2. This solid line clearly shows 
the expected behavior of f{x): f{x) = 1 for a; ^ 0, 
and f{x — > oo) cx x, with a crossover in between. The 
crossover regime can be estimated by setting B = A. The 
relaxation time can be inferred from the ab plane resis- 
tivity measurement iii at B=10 T, p ^ 20/ifi • cm, which 
gives T ^ 10~^^ sec. So, for A = O.leV, corresponding 
to Xc — X ~ 0.005 if we take A(x) = 0.7 y/l — x/xcbV, 
the crossover is B ^ ^qHA/ (vpr) ^ lOOT, where 
vp 4eV-A. For these parameters, the breakdown field 
is slightly larger. 



VI. SENSITIVITY TO PARAMETERS 

In this section we examine the sensitivity of our results 
to different choices of parameters. 

Let us first consider the x < Xc region. FigElshows nu- 
merical results of Rh for parameter choices: ti = 0.38eV, 
t2 = 0.32ii, ^3 = 0.075eV and A{x = 0) = 0.6eV (dashed 
line), and = 0.32ti, = 0.5^2 and A{x = 0) = 0.5eV 
(solid line), which can be compared with the solid and 
dashed lines in Fig|2| All of these curves share the same 
feature: the existence of xi. In FigO xi « 0.12, while in 
Figini xi « 0.09 (dashed line) and xi ~ 0.06 (solid line). 
For X > xi, both electron pocket (centered at (0,7r)) and 
hole pockets (centered at (7r/2, ±7r/2)) are present, and 
as doping decreases, both of them contribute to the de- 
crease of Rh, while for x < xi, the hole pockets vanish 
(the E~ band is filled), and thus only electron pocket 
contributes to Rh, which leads to the slowing down of 
decrease in Rh as decreasing doping. For x < xi, the be- 
havior of Rh is quite universal, approximately ~ — 1/x 
(the Hall coefficient for free electrons, with x electron 
density), as shown in Fig 13 Different choices of param- 
eters change the actual value of xi. The numerical am- 
plitude of Rh{xi) is approximately l/xi, by continuity. 
There is still another possibility that if the opening of 
SDW gap is not fast enough, the hole pockets survive for 
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FIG. 6: Numerical results of Hall coefficient in SDW ordered 
state for different choices of parameters. Dashed line: t\ = 
0.38eV, t2 = 0.32ii, is = 0.075ey and A(a; = 0) = G.&eV. 
Solid line: ts = 0.5t2, t2 = 0.32ii, A(a:: = 0) = 0.5eV. Inset: 
Rh for A{x = 0) = O.SeF, with the same band parameters as 
the solid line in the main panel. 



X down to 0. In this case, there is an upturn when x 
is smaller than some value, with a much smaller relative 
change in the magnitude of Rr, as shown in the inset of 
Fig El We do not believe this possibility is relevant to 
the electron doped cuprates. 

In the regime x > Xc, we can change band parameters 
t2 and ^3 in Eq , and calculate Rh for these choices of 
parameters. Fig [7| shows the results for different choices 
of band parameters: ^2 — 0.32ti and ^3 — 0.075eV (heavy 
solid line); ^2 = 0.2ii and ^3 = 0.5^2 (dashed line); 
t2 = 0.45^1 and ^3 = 0.5^2 (dotted hue). The Fermi 
surfaces corresponding to these choices of parameters at 
X = 0.16 are shown in Fig|Sl It can be seen that (i) dif- 
ferent choices of band parameters do not lead to drastic 
difference in Rh] (ii) for each choice of band parameters, 
Rh is quite insensitive to change in doping, in contrast 
to the experimental data, shown as solid circles; (iii) our 
numerical Rh are close to 1/(1 — x) (the Hall coefficient 
for circular Fermi surface, with 1 — x hole density), shown 
as light solid line in Fig \J\ Both our numerical results 
and l/j^l — x) are not close to the experimental data of 
Ref. 11. The difference is largest at x = 0.19, the point 
farthest from criticality, where one would expect the dif- 
ference to be smallest. 



VII. X > Xc 

Fig H shows that the measured Rh at x > Xc shows 
a pronounced doping dependence which is incompatible 
with the mean field calculation at fixed band parameters. 
Fig 13 also shows that a variation of the correct order of 
magnitude can be obtained if a doping dependence of the 
band parameters is assumed. In this section we investi- 
gate whether this doping dependence can be viewed as a 
precursor of the spin density wave instability by calculat- 
ing the effect of spin fiuctuations on the Fermi surface. 
This is given by the real part of the self energy. The 



FIG. 7: Numerical results of Hall coefficient in normal state 
for different choices of band parameters. Heavy solid line: 
t2 — 0.32fi and t-j, — 0.075eF; dashed line: t2 — 0.2ti and 
ts — 0.5t2; dotted line: t2 — 0.45fi and ts = 0.5t2; light solid 
line: curve 1/(1 — x); solid circles: experimental dataJ^ 



X 




r 



FIG. 8: Normal state Fermi surfaces in the first quadrant 
of Brillouin zone for different choices of band parameters at 
X — 0.16. Solid line: t2 — 0.32ti and ts — 0.0756^; dashed 
line: t2 — 0.2ii and ts = 0.5t2; dotted line: £2 = 0.45ii and 
ta = 0.5t2. 

leading approximation to this is: 

S(k, iLu) = -g^T J (dq) ^ ^(k + q,iLU + zr!„)I?(q, ifln) 

(34) 

where Q and V are electron and spin fluctuation Matsub- 
ara Green's functions, respectively. In our calculation, 
we shall take: 0{p,iuj) — {iuj — ep)~^ and I?(q, zi7„) — 
-(r, + |r!„|)-i, with r, = r(r + ^2(q_Q)2)^ whereris 
the energy scale characteristic of spin fiuctuations, ^0 is 
coherent length, r characterizes the distance to critical- 
ity, and Q = (tt, tt). 

This self energy for icu — can be integrated immedi- 
ately, at T ^ 0, if both ek+g/Eo {Eq is an energy cutoff) 
and r are small. The result may be represented as 

S(k, ^u^0)^ -0.5Ae.+g In . 2 + ( ' /j;„)4) (35) 
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where, A is a dimensionless constant. 

We can take Sp + E(p,za' = 0) as our new electronic 
energy dispersion relation, and integrate Eqs |14I16| nu- 
merically. The Hall coefficient Rh as function of r, for 
A — 0.2 (solid line) and A — 0.4 (dashed line), is shown in 
Fig 1^1 where the parameters are taken as: ti = 0.38eV, 
t2 = 0.5ti, = 0.075eV, and Eq = 0.88eV. Because 
the relationship between r and doping has not yet estab- 
lished, in our numerical calculation, we keep the chemical 
potential fixed at the unrenormalized value for x = 0.16. 
FigllOlshows Fermi surfaces for the case A = 0.2, r — 0.01 
and the case A = 0.4, r = 0.1, as well as the original un- 
renormalized one. Both of these Fermi surfaces give Rh 
about 25% less than the unrenormalized value. 




FIG. 9: Numerical result for Hall coefficient as a function of 
r when electronic energy dispersion is Sp + E(p,iaj = 0), for 
coupling constant A — 0.2 (solid line) and A — 0.4 (dashed 
line). Other parameters are: t\ = 0.38el/, t2 = 0.5ti, t^ — 
0.075eF, the energy cutoff Eo = 0.88eV. 
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FIG. 10: Fermi surfaces in normal state, with renormalized 
energy dispersion relation, Eq \!I5( : A = (solid line); A = 0.2 
and r = 0.01 (dashed line); X — 0.4 and r — 0.1 (dotted line). 
The numerical Rh corresponding to the dashed line is about 
20% less than that corresponding to the solid line, and Rh 
corresponding to dotted line is about 30% less than that of the 
solid line. 



VIII. CONCLUSION: DISCUSSION OF 
EXPERIMENT 

We have used mean field theory to calculate the 
changes in the Hall resistivity of a two dimensional sys- 
tem which undergoes a spin density wave instability. 
The model studied was chosen to be appropriate to the 
electron-doped cuprates. The important physics is that 
as X is decreased below a critical value Xc, the Fermi 
surface is rearranged (Fig^, leading to changes in Rh- 
These changes were discussed, in the context of a theory 
of Cr, by Bazaliy et al,^^ and the related discussions, 
in the context of the d-density wave model of cuprates, 
were given by Ref. 9. The new results presented here 
include consideration of a wide range of dopings and gap 
values (FigO, a systematic treatment of all relevant ef- 
fects (Sec. ini|l . a discussion of the weak field-strong 
field crossover (Sec. IIV|l . a treatment of precursor ef- 
fects (Sec. IVn|l . and detailed application of the model to 
the electron-doped cuprates (Fig 13 and conclusion). We 
note, unfortunately, that the specifics of the critical be- 
havior and weak to strong field crossover will be visible 
(at presently accessible fields B < SOT) only at samples 
doped very close {Sx ~ 0.005) of the critical doping, and 
at the present level of sample purity, these effects will be 
complicated by magnetic breakdown. The data of Ref. 
im suggest Xc ~ 0.165 while the data of Ref. indicate 
xi > 0.1 and is most probably ~ 0.12. Thus the gap ap- 
parently opens quite fast with doping. However, impor- 
tant discrepancies remain between data and calculation. 
At X < xi our calculation shows Rh ~ 1/a; (as expected 
for a Fermi surface with one electron pocket containing 
X electrons). The data of Ref. ^12, are qualitatively con- 
sistent with this trend, but are offset, suggesting that 
the actual doping is slightly lower than the nominal one. 
The X = 0.11 data point of Ref. JT, is qualitatively in- 
consistent with theory. Also, the magnitude reportedii 
of Rh at X > Xc is inconsistent with theory. One expects 
that as X is increased beyond Xc, Rh should revert to the 
band value, Rh ^ 1/(1 — x). While the sign and order of 
magnitude are correct, the numerical value is too large; 
there is no reasonable choice of band parameters which 
reproduces the magnitude to better than a factor of two. 

Our work suggests several important directions for fu- 
ture research. On the experimental side, we suggest that 
our finding that Rh ~ 1/x for x < xi can be used to 
determine the actual doping of lightly doped samples 
{x < 0.1). We also note that experiments which pinned 
down the behavior in the Xc > x > xi regime and more 
precisely located xi would be very helpful. Finally, clar- 
ifying the behavior at x > Xc is potentially of consid- 
erable interest, as it gives (in principle) evidence of the 
modification of the Fermi suface by finite range SDW cor- 
relations. On the theoretical side, important directions 
include combining the present considerations with those 
of Ref. 10 to obtain a theory of the Nernst effect, and im- 
proving the theoretical treatment of the SDW precursor 
effects. Work on the AC Hall effect is in progress. 



9 



Acknowledgements We thank R. Greene, Y. Dagan and and the possibiHty of magnetic breakdown. This work 
M. R. Norman for helpful discussions and A. J. Schofield was supported by NSF-DMR-0431350. 
for drawing our attention to the strong-field crossover 



^ S. Sachdev, Rev. Mod. Phys. 75, 913 (2003). 

^ P. Monthoux, A. V. Balatsky, and D. Pines, Phys. Rev. 

Lett. 67, 3448 (1991). 
^ S. A. Kivelson, I. P. Bindloss, E. Fradkin, V. Oganesyan, J. 

M. Tranquada, A. Kapitulnik, and C. Howald, Rev. Mod. 

Phys. 75, 1201 (2003). 

C. M. Varma, Phys. Rev. Lett. 83, 3538 (1999). 
^ G. M. Luke et. al., Phys. Rev. B 42, 7981 (1990). 
® M. Matsuda, Y. Endoh, K. Yamada, H. Kojima, L Tanaka, 

R. J. Birgeneau, M. A. Kastner, and G. Shirane, Phys. 

Rev. B 45, 12548 (1992). 

P. K. Mang, O. P. Vajk, A. Arvanitaki, J. W. Lynn, and 
M. Greven, Phys. Rev. Lett. 93, 027002 (2004). 

* M. R. Norman, Q. Si, Ya. B. Bazaliy, and R. Ramazashvili, 
Phys. Rev. Lett. 90, 116601 (2003). 

® S. Chakravarty, C. Nayak, S. Tewari, and X. Yang, Phys. 



Rev. Lett. 89, 277003 (2002). 
° V. Oganesyan and I. Ussishkin, Phys. Rev. B 70, 054503 
(2004). 

1 Y. Dagan, M. M. Qazilbash, C. P. Hill, V. N. Kulkarni, 
and R. L. Greene, Phys. Rev. Lett. 92, 167001 (2004). 

^ Y. Onose, Y. Taguchi, K. Ishizaka, and Y. Tokura, Phys. 
Rev. Lett. 87, 217001/1-4 (2001). 

^ R. G. Chambers, Proc. Roy. Soc. (London) A238, 344 
(1957). 

E. L Blount, Phys. Rev. 126, 1636 (1962). 

5 N. P. Ong, Phys. Rev. B 43, 193 (1991). 

® Ya. B. Bazaliy, R. Ramazashvili, Q. Si, and M. R. Norman, 

Phys. Rev. B 69, 144423 (2004).^ 

J. Fenton and A. J. Schofield, co nd-mat/0507245| 



